Stability in the energy space for chains of solitons of the 
one-dimensional Gross-Pitaevskii equation 
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Abstract 



We establish the stability in the energy space for sums of solitons of the one-dimensional 
Gross-Pitaevskii equation when their speeds are mutually distinct and distinct from zero, 
and when the solitons are initially well-separated and spatially ordered according to their 
speeds. 

1 Introduction 

' In this paper, we focus on the one-dimensional Gross-Pitaevskii equation 

*S '• id t t> + 8%® + V(l-\t>\ 2 ) =0, (GP) 

for a function $ : 1 x I 4 C, supplemented with the boundary condition at infinity 

\^(x,t)\ — > 1, as |x| — > +00. (1) 

The Gross-Pitaevskii equation is a defocusing nonlinear Schrodinger equation. In one spatial 
dimension, solutions with finite Ginzburg-Landau energy 

o 



U 1**1" + i /.»-!• 



2\2 



are globally defined, and the energy is conserved along the flow. 

In 1972, V.E. Zakharov and A.B. Shabat |20| have shown that the Gross-Pitaevskii equation 
supplemented with condition (JTJ) is integrable by means of an inverse scattering transform. As 
a consequence, they exhibited exact soliton and multi-soliton solutions. It is generally believed 
that such solutions play a central role in the long-time asymptotics of the flow, and sometimes 
conjectured that any initial datum with finite Ginzburg-Landau energy should eventually resolve 
into a multi-soliton solution plus a vanishing (in an appropriate sense) dispersive tail. Detailed 
asymptotic expansions for smooth and fast decaying solutions, based on their scattering data 
at the initial time, have actually been obtained by A.H. Vartanian in |19j . using reductions to 
Riemann-Hilbert problems. To our knowledge, for initial data in the energy space, there is no 
rigorous evidence supporting the aforementioned behaviour further than the orbital stability of 
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a single soliton (see e.g. |12| [TJ [2j [TO]). The main goal of this paper is to provide a partial 
justification through the proof of the orbital stability in the energy space of some finite sums of 
solitons along the Gross-Pitaevskii flow. Our approach does not make any use of the integrability 
of the equation, and is largely influenced by the corresponding study of the generalized Korteweg- 
de Vries equation (see |15| ) and of the nonlinear Schrodinger equation with vanishing data at 
infinity (see [16] ) by Y. Martel, F. Merle and T.-P. Tsai. The main ingredient is a monotonicity 
formula for the momentum (see Subsection I3.2|) . 



1.1 Soliton solutions to the Gross-Pitaevskii equation 

A soliton with speed c is a solution to (|GPp which takes the form 

*(x,t) := U c (x - ct). 
Its profile U c is a solution to the ordinary differential equation 

- icd x u c + d 2 xx u c + u c (i - \u c \ 2 ) = o. (2) 

The solutions to ([2]) with finite energy are classified and indeed explicit. For |c| > \/2, all of them 
are identically constant. When |c| < \/2, there exist non-constant solutions with finite energy. 
Up to the invariances of the problem, i.e. the multiplication by a constant of modulus one and 
the invariance by translation, they are uniquely given by the expression 



<•■"■.)• \ . c 



Uc(x) = \l =-^- ( (3) 

Notice that the soliton U c does not vanish on R when c / 0. As a matter of fact, our stability 
analysis below requires non- vanishing solutions. Indeed, it is performed on a reformulation of 
([GPP which only makes sense for such solutions, and which we introduce next. 

1.2 Hydrodynamical form of the Gross-Pitaevskii equation 

Provided a solution \& to (|GP|) does not vanish, it may be written, at least formally, as 

^ := gexpiip, 

where g := |^|. Still on a formal level, the functions r\ := 1 — g 2 and v := d x (f are then solutions 
to an hydrodynamical form of (|GPp . namely 

d t r\ = d x (2v - 2rjv) , 



.2(1-7/)/ 4(1 -T]f 

For k > 0, we introduce the Hilbert spaces X k (R) := H k+1 (R) x H k (R), equipped with the norm 



:— ll 7 ?ll/i'' s + 1 \\ v \\H k > 



and their open subsets 

MV k (R) := \(r},v) G X k (R), s.t. maxrj(x) < l), 



1 The notation J\fV stands for non- vanishing. 
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consisting of states with non-vanishing corresponding ^. In the sequel, we set X(R) := X (R) 
and MV (R) :=AfV°{R). 

The Ginzburg-Landau energy 

E(r,,v):= [ e(r,,v):=lf ^ + \ I (1 - ^ + \ [ V \ (4) 
Jr 8 J R 1 - n 2 J R 4 J m 

and the momentum 

P(r),v) := / p(?7,u) := - / rjv, 
JR z JR 

are well-defined and smooth functionals on A/"V(M). In view of (£[]), the space AAV(M) is the 
energy space for states (r],v) satisfying the non-vanishing condition max xe jR r](x) < 1. 

Concerning the Cauchy theory for (jHGPp . we have 

Theorem 1 ([18J). Let k G N and (n°,v ) G MV h (R). There exists a maximal time T max > 
and a unique solution (r],v) G C°([0, T max ), AAV fc (M)) to equation ()HGP|) with initial datum 
(r]°,v°). The maximal time T max is continuous with respect to the initial datum (t]°,v ), and is 
characterized by 

lim maxnfx) = 1 if T max < +oo. 

t^T max x£R 

The flow map (n ,v ) i— > {i],v) is locally well-defined and continuous from AfV k (M) to C°([0,T], 
MV k (R)) for any T < T max . Moreover, the energy E and the momentum P are constant along 
the flow. 

We refer to |18j for a proof. It is a consequence of the Cauchy theory for (|GPp in Zhidkov 
spaces Z k {M) (see e.g. |21l 18) [9]). and the fact that the mapping ^ i— > (n,v) is locally Lipschitz 
from Z k (R) to X k (R) provided that (rj,v) G MV(R). The uniqueness part requires slightly more 
care. 



1.3 Statement of the main result 

When c ^ 0, the solitons U c described in Subsection 11.11 do not vanish and may thus be written 
under the form 

U c := g c expiip c . 

In view of formula ([3]), the maps r\ c := 1 — and v c := d x (p c are given byH 

/ \ 2-c 2 cri c (x) c(2 - c 2 ) 

^cl^j = ; ? — ^ and vAx) = -. — T = ; — k . (5) 

2ch (^ra ) 2 2(1 - 77 C (X)) 2 ( 2c h(^a) 2 _ 2 + c 2) 



2 When may be written as ^ = gexpiip, the quantities £ (^) and v), with 77 = 1 — g 2 and n = d x tp, are 
exactly equal. 

3 A striking point in formula ([5jl is the fact that, up to some scaling, the map r/ c is exactly the soliton of the 
Korteweg-de Vries equation 

d t u + dlu + ud x u = 0. (KdV) 

Like the Gross-Pitaevskii equation, the Korteweg-de Vries equation indeed owns a family of solitons u a , given, up 
to translations, by 

u a (x) ■— — , 

for some positive speed a. Actually, the link between the two equations does not reduce to the level of the solitons. 
Both of them are integrable by means of the inverse scattering transform (see e.g. [17] for the Korteweg-de Vries 
equation), and moreover, as was pointed out in 2, 3 (see also [5]), the Korteweg-de Vries equation appears as 
the limit of system (|HGP[) in some long-wave regime. 
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In the sequel, we set 

Qc,a ■= {ilea, V c ,a) ■= {Vc(- ~ o), V c (- - a)) , 

for < |c| < y/2 and a£l. We introduce the set of admissible speeds as 

AdniAr := {c := (a,...,c N ) G (~V2,V2) N , s.t. c k / 0, VI < k < N}, 

and, for c G Adm^v and o := (a\, . . . , aAr) G R , we define 

N , N N \ / ^ N \ 

Rc,a-= {Vc,a,Vc, a ) = ^Qc k ,a k = ( ^2Vc k ,a k ,^2v Cktak ) = ( 77 Cfc (• - Oft), U Cfc (• - Ofc) J . 

fe=l ^Jfc=l fc=l ' ^fc=l fc=l ' 

In the sequel, we refer to small perturbations of the functions i? ca as chains of solitons. We also 
denote 

v< : = minjy^ - c|, 1 < k < n\. 

Our main result is 

Theorem 2. Let c* G Admjy be such that 

cl<...<c* N . (6) 

There exist positive numbers a*, L* and A*, depending^ only on c* , such that the following 
holds. Assume that (r] ,v ) G J\fV(M) satisfies 



a :-- 



(r,°,v°)-R c * a0 



x 



< a 



for some points a = (a®, . . . , a%) G M N such that 

L° := min {a° k+1 - a° k , 1 < k < N - l} > L*. 



(7) 



(8) 



Then the unique solution (77, u) to (jHGPj) with initial datum (r/°, u°) is globally defined on [0, +00) 
and f/iere exists a function a = (ai, . . . , ajv) G C 1 ([0, +00), WL N ) such that 



N 



£>'*(*) -4| < A*(a°+exp 



fe=l 



33 



and 



fa(-,t),v(-,t))-flc*,a 



(t) 



< ^*(a° + exp 



x 



v t *L> 
33 



(9) 



(10) 



/or any i G [0, +00). 



In different terms, chains of solitons with strictly ordered speeds are preserved along the 
Gross-Pitaevskii flow in positive time. Since the initial speeds are strictly ordered, the mutual 
distances between the individual solitons appearing in the chain increase proportionally with 
time, so that the interactions between solitons decrease to zero. As a matter of fact, integrating 
© yields 

a k+1 (t) - a k (t) > (c* k+1 -c* k - 2 A* (a + exp ( - ^~)))t + a*+i(0) - a fc (0). 

4 A close inspection of the proof actually shows that a*, L* and A* only depend on a lower bound on v c * , on 
min{|c* — ell, 1 < j < k < N} and on min{|c£|, 1 < k < N}. 
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A typical example of chain of solitons is given by the exact multi-solitons constructed by Zakharov 
and Shabat |20] (see also [7J). In particular, Theorem [2] provides the orbital stability of multi- 
solitons in the energy space in positive time. 

In view of the time reversibility of the Gross-Pitaevskii equation, Theorem [2] also provides the 
stability in negative time of chains of solitons with reversely ordered speeds. To our knowledge, 
the question of the stability in the energy space of chains of solitons for both positive and 
negative times remains open when the chain includes at least two solitons. For the generalized 
Korteweg-de Vries equations, results of this type have recently been obtained by Martel and 
Merle [131 H3]. 

Finally, we believe that our arguments still apply to non-integrable versions of the nonlinear 
Schrodinger equation with non-vanishing condition at infinity, provided individual solitons are 
themselves orbitally stable. The latter condition has recently been treated in details in a work 
by Chiron [3], where numerous examples are presented. 

1.4 Elements in the proofs 

We will present a detailed description of the elements in the proof of Theorem [2] in Subsections 
[A] to [G] below. We precede this description by a voluntarily oversimplified overview in the case 
of a two-soliton, in the hope that this may ease the reading of the detailed description. 

The solitons Q c may be characterized as minimizers of some variational problems. In terms 
of Q c = (rj c ,v c ), equation ([2]) reduces to 

E'(Q c )+cP'{Q c ) = 0. (11) 

Equation (|lip is the Euler-Lagrange equation for the minimization of the Ginzburg-Landau en- 
ergy E under fixed scalar momentum P, and the speed c is the corresponding Lagrange multiplier. 
As a matter of fact, we have 

E(Q C ) = mm {e(t],v), (rj, v) G X(R) s.t. rj < 1 and P(rj, v) = P(Q C )}, (12) 

for any c G (— \/2, V%) \ {0} (see e.g. |lj), and Q c is the only minimizer up to the invariances 
of the problem. This minimizing property, combined with the conservation of E and P by the 
Gross-Pitaevskii flow, yields the orbital stability of single solitons. Instability of a chain of 
solitons, should it happen, could thus only arise from their mutual interactions. Assume for 
simplicity that at some fixed time the solution at hand is a chain of two solitons 

(r),v) — Qci,ai 4" Qc2,02) 

where a 2 — d\ S> 1. In particular, the total energy E and momentum P satisfy 

E~E(Q C1 )+E(Q C2 ) and P ~ P(Q Cl ) + P(Q C2 ). 
Assume next that at some later time, interaction has occurred and 

(r l ,v)~B l (--b 1 ) + B 2 (--b 2 ), 
where still b 2 — b\ 3> 1, but B\ and B 2 need no longer be solitons. Set 

6P:=P(Q C1 )-P(B 1 ). 
Since the total momentum is an invariant of the flow, we also have 

5P^P(B 2 )-P(Q C2 ), 
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so that 5P may be understood as the amount of momentum transfered from Q C1 ai towards Q C2 ,a 2 
(we let aside here the possibility of transfer to the background) . 

Prom (|12p . we infer that the energies E(B\) and E(B 2 ) are necessarily larger than the energies 
of the corresponding solitons with momentum P{B\), respectively P(B 2 ). Since, in particular 
from (HH), 



we obtain at first order in 8P, 
and 



dc 

E(Bi) > E(Q C1 ) + c x 8P (13) 



E{B 2 ) > E(Q C2 ) - c 2 5P. (14) 
On the other hand, by conservation of the total energy, we have 

E(Q C1 ) + E{Q C2 ) ~ E{B{) + E(B 2 ), 

so that, by (THj) . 

E(B X ) < E(Q C1 ) + E(Q C2 ) - E{B 2 ) < E(Q Cl ) + c 2 5P. (15) 
Combining (I13p and ()15p yields 



(c 2 - ciJJP > 0. (16) 

The last and key observation is a monotonicity formula, which states that in the above configu- 
ration necessarily 

6P < 0. (17) 

Therefore, if the speeds are ordered as in Theorem [21 namely c 2 > c±, then (fl~6|) may hold only 
if 5P = 0. Orbital stability of the chain then follows from the orbital stability of the individual 
solitons. 

Remark 1. The monotonicity formula (fT7|) has a simple physical interpretation. Momentum is 
a signed quantity which, in the present normalization, is non-negative for "waves" or "particles" 
travelling to the left and non-positive for those travelling to the right. Since Q C1 ai is located at 
the left of Qc- 2 ,a 2 (recall that a 2 —a\ 3> 1 by assumption), the momentum transferred from Q C i,ai 
towards Q C2 .a 2 (i- e - needs to travel from the left to the right, and is therefore non-positive. 
Equivalently (and alternatively from the point of view of Q C2 ,a 2 ), the momentum transferred 
from Q C2 ,a 2 towards Q Cli ai) i- e - — SP, needs to travel from the right to the left, and is therefore 
non-negative. 

We are now in position to present the detailed description of the elements in the proof. For 
later reference, we mention the explicit values 

P(Q C ) = arctan f— L= ) + - — - sign(c) J, (18) 



and 



E{Qc) = $-4^' (19) 



To be more precise, the transferred momentum needs to travel from the left to the right in a reference frame 
attached to Q ci , ai , i.e. moving at speed c\. Dispersive waves have a speed larger or equal in absolute value than 
the speed of sound y/2 given by the dispersion relation. Since solitons are subsonic, a dispersive (hence sonic 
or supersonic) wave going to the right in the reference frame attached to a soliton also goes to the right in the 
frame of the lab. This peculiarity is a major difference between the Gross-Pitaevskii equation and the nonlinear 
Schrodinger equations with zero density at spatial infinity, where the speeds of solitons and dispersive waves may 
overlap. 



A Minimizing properties of solitons 

We will rely on a quantitative version of the minimizing property (|12p . Define the quadratic 
form 

H c :=E"(Q c )+cP"{Q c ). 

In explicit form, 

where we have written £ = (e„,e v ). We have 
Proposition 1. There exists A c > such that 

H c (e) > A c \\e\\x, (21) 

for any e € X(R) satisfying the orthogonality conditions 

(e,d x Q c ) L 2 m =P'(Q c )(e) = 0. (22) 

Moreover, A c is uniformly bounded from below by a positive constant for c in a compact subset 
o/(-V2,V2)\{0}. 

Remark 2. The first orthogonality condition in (|22l) is related to the invariance with respect to 
translations of (|HGP]> . It may be rephrased according to the equality 

{e,d x Q c ) L 2 = 1 — ^|| - Q c (- = 0. 

It is possible to choose the condition involving the momentum P as an orthogonality condition 
due to the property 

|(P(Q c ))/0, 

which follows from the explicit value (|18|) of P(Q C ). 

In particular, if £ := (rj,v) — Q c satisfies the orthogonality conditions ([22|) . then 

(E + cP)( V ,v) > (£: + cP)(Q c ) + A c ||£|^ + C»(||£|^), (23) 

as ||e||x - * 0. 

Our next goal is to obtain a coercivity inequality in the spirit of (|23p for functions (r/, v) which 
do not belong to a small neighbourhood of a single soliton but instead to a neighbourhood of a 
finite sum of well-separated solitons like the one described by j7|) and ([8]). We begin with the 
treatment of the orthogonality conditions in that case. 

B Orthogonal decomposition of a chain of solitons 

Given a positive parameter L, we define the set of admissible positions as 

Posat(L) := {o := (ai,.. . ,a N ) G R N , s.t. a k+1 > a k + L, VI < k < N - l}, 
and, for a > and c = (ci, . . . , c^) £ Adm^r, we set 

U c (a,L) := \(r),v) e X(R), s.t. inf II (77,1;) - a \\ < a\, 
I oePos]v(L) ' J 
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and 

/j c := min {|cfc|, 1 < k < A^}. 

When a is small enough and L is sufficiently large, we claim that the set U c (a,L) provides a 
suitable framework to perform satisfactory decompositions. More precisely, we have 

Proposition 2. Let c* G Adm^r. There exist positive numbers ct\, L\ and K\, depending 
only on c* , and two continuously differentiable functions £ G C 1 ^* (ai, L\), Adm^r) and 21 G 
C 1 (Z^ C * (ai, Li), R ) swc/i that for any (77, u) ^ U c *(a\, L\), the function 

e ■- ( v ,v) - R c , a , 

where c = (ci,...,Cjy) := £(17, v) and a = (ai, . . . , Ojv) := 2l(r/, v), satisfies the orthogonality 
conditions 

(e,d x Q Ck , ak ) L 2 = P'{Q Ck>ak ){e) = 0, (24) 
/or any 1 < A; < AT. Moreover, if 

\\{r),v) - Rc*, a *\\ x < a, 
for some a* = (a*, . . . , a^) G Posat(-^)> mtt L> L\ and a < a\, then 

N N 

\\ £ \\x + ^\ck- 41 + ' afe ~ - ^ 1<x ( 25 ) 

fc=l fc=l 

Increasing Li and decreasing a,\ in the statement of Proposition [2] if necessary, one may further 
infer from (|25p the following facts which we will use throughout. 



Corollary 1. Under the assumptions of Proposition^ whenever (n,v) G V( c *(ai,L) with L > L\, 
it follows that (rj,v) G NV(M), since 

1-V>^- (26) 

Moreover, we have 

2t(77, v) G Posat (L — l) , (27) 
A%M) > ■ (29) 



C Almost minimizing properties close to a sum of solitons 

Given c* := (c*, . . . ,c* N ) G Admjv and (r),v) G U c *(ai, L\), we write 

(r?,u) = R c , a + e, 

according to Proposition [2j where e satisfies the orthogonality conditions (|24p . Our next goal 
is to show that the sum i? ca possesses almost minimizing properties. Since each soliton Q Ck ,a k 
has its own speed c/%, the function E + cP, whose coercivity properties were exhibited in ([23|) in 
the case of a single soliton Q c , has to be replaced by some function which, close to each location 
point a/-, resembles E + c^P . 
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More precisely, we define the functions [^| 

' 1 







We next set 



where 



o-k-1 + a k 



N 



if fe = 1, 

if 2 < k < N, 
if jfe = JV + l. 



F(r/, u) := E(tj, v) + ^ c fe P fc (??, u), 



fe=i 



/ p(T7,tO(¥fc-tffc+i) 



Notice that by construction, 



N 



(** - tffc+i) = *i - = i, 



*!=i 

and that, provided the points a k are sufficiently far apart one from each other, we have for each 
k, — — 1 in a neighborhood of Q Ck ,a k and ty k — ^>k+i — in the neighborhood of Q Cj , aj , 

for j ^ k. Notice also that, even though it does not explicitly appear in the notation, ty k depends 
on (rj, v) through the points a k . This will be of particular importance when we introduce time 
dependence later on. 

In order to estimate the function F, we localize the function e according to the following 
decomposition. Let < r < z/ c */16 to be defined later (see (|38|) ). depending only on c*. We set 



$k(x) 



1 



for 1< k < N, and 



k,k+l 



- (It (x-a k + — 



L1 

t \ x - a k - — 



Ul+(T(x-a N -^ 



(t(x - a k+1 + If 



if k = 0, 

if Kk< N, 

if k = N. 



Once more, we note that by construction, 

N 



N 



k=l k=0 

and that each term in the sum is non-negative. We finally define 



for 1 < k < N, and 
for < k < N. 



e k := $ fc (- + a fc )2e(- + a k ), 
i 



(30) 

(31) 
(32) 



6 In the definition of one could a priori replace the constant z/ c »/16 by an arbitrary positive constant v* . 
In order that the monotonicity formula later stated in Proposition [5] holds, the choice of v* is actually restricted. 
Up to a constant multiplicity factor, ^ c */16 is presumably optimal. 
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Proposition 3. Assume that (rj,v) = i? Ci(J + e G L( c *(ai, L), with L > L\. Then, 

N N N 

E(r,,v) =J2 E (Qc k ) + -^E"(Q Ck )(e k ,e k ) + J^>"(0)( £fc 

fc=i fe=i k=o (33) 

+ o( (r + exp ( - t%)) \\s\\ x ) + o( \\e\\ x ) + ofiexp / "' ' 



2 

and 

P k (v,v) =P{Q Ck ) + \p"{Q Ck ){ek,e k ) + -i^'(0)(e fc , fc+ i,e fc , fc+ i) + ^i^'(0)(e fc _i >fcj e fc _ lifc ) 

2 (34) 



+ o(M _^ ))+o(exp( _^ )lk 



Remark 3. Here and in the sequel, we have found convenient to simplify the presentation of 
our estimates by introducing the notation O. We adopt as a definition that one is allowed to 
substitute a quantity S by the notation 0{T) if and only if there exists a positive constant K, 
depending only on c* , such that 

\S\ < KT. 

Later on, we will analyze the evolution in time of the evaluation of F along a solution (77, v) 
of ()HGP[) . The analysis is complicated by the fact that the speeds c k are themselves functions 
of (ji,v). To circumvent this difficulty, it turns out that it is sufficient to consider instead of F 
the function G defined on U c * {a\,L{) by 

N 



k=l 

Notice here that the speeds are the reference speeds ct. Concerning the functional G, we have 
the following consequence of Proposition [3J 

Corollary 2. Assume that (rj,v) = R c ^ a + e £ U c *(ai,L), with L > L\. Then, 

N N N N 

G( V ,v) = £ ( E (Qc%) + 4P(Q C% )) + -Y,H Ck (e k ) + -Y, H o(zk, k+ i) + ^(E ^ " c l 

fc=l fc=l k=0 k=l 

+ C((r + exp ( - ^)) \\e\\ x ) +o(\\e\\ x ) + o(jDe*p ( - ^ 
where we have set 

#§(eo,i) = ( J B ,/ (0) + ciP[ / (0))(eo ) i,eo,i) ) 

H*(e kik+1 ) = {E"(0) + c k P',!(0) + c k+1 P!! +1 (0))(e k:k+1 ,e ktk+1 ), VI < fc < JV - 1, 
Hq ( £ n,n+i) = (E"(0) + cnP'n(0))(£n,n+i,£n,n+i)- 

It is not possible to invoke Proposition[T]directly to obtain a lower bound on the terms H Ck {e k ) 
since only satisfies the orthogonality conditions (|22[) asymptotically as Li — )• +00. We obtain 



Lemma 1. Under the assumptions of Corollary^ there exists a positive constant A* , depending 
only on c* , such that 

H Ck (e k ) > A* ||e fc f x + O (l\ exp ( - ^) \\e\\ x ) , (35) 
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for any 1 < k < N, and such that 

^o( £ fc,A:+l) > A-*[| £ fc,fe+l||x! (36) 

/or any < k < N . Moreover, 

N N 

^ \\ £ k\\x + 5Z - ll e llx- ( 37 ) 

fc=l fc=0 

In view of Corollary [2] and Lemma [H we now fix the value of r. sufficiently small, and then 
increase the value of L\ , if necessary, in such a way that 

0((r + Lle, P (- T -^))\\ef x )<^\\sf x . (38) 
With these choices we obtain 

Corollary 3. Assume that (n,v) = i? Ci0 + e G U c *(a\, L), with L > L\. Then, 

N A * N 

G( V , v)>Y, {E(Qci) + 4P(Q C *)) + — \\e\\ x + 0(\\e\\ x ) + o( £ \c k - c* k \ 2 ) 

fc=l fe=i (39) 

as well as 

N N 

G( V ,v) <^( J E(Q c .) + 4P(Q c *))+0(|| e |||)+0(^|c fc -4| 2 )+0(Lexp(-^ F )). (40) 

k=l ' fc=l 

Notice that, up to this point, the entire analysis has been independent of time. In the next 
section, we combine the previous results with information extracted from equation (jHCPp . 

D Dynamics of the modulation parameters 

In this subsection, as well as in Subsections [El and IF1 below, c* G AdrriTv is fixed and we assume 
that T > 0, and (r),v) G C°([0,T],X(M)) is a given solution to equation (jHGPj) such that for 
some a < a\ and L > L±, 

{r,(;t),v{;t)) £U c *(a,L), 

for any t G [0,T]. 

For t G [0, T], we define 

c(i) = (ci(i), . . . , civ(t)) := £(«(-, t), 7/(-, t)) and a(t) = (oi(t), . . . , a N (t)) := %{v(;t), V (;t)) , 
and 

£ 0)*) : = {v(-,t),v(-,t)) -R c ( t ), a (t), 

where the functions £ and 21 are given by Proposition [2 

The following proposition expresses the fact that the modulation parameters should follow 
those of the underlying solitons, as long as the solution remains close to them. 



7 By the left-hand side of (|38[) . we actually mean all the terms which we have previously estimated in Corollary 
[2] and Lemma [T] as the corresponding O. 
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Proposition 4. There exist positive numbers a% < ol\ and L2 > L±, depending only on c* , such 
that, if a < (X2 and L > L2, then c G (^([O, T], Adm/v), a G C 1 ([0, T], 1^), and we have the 
estimates 

N 

£ (|a' fc (i) - c fc (t)| + = 0(lk(-,t)||x) + ©(Lexp ( - ^-)), (41) 

fc=i 

/or any t G [0, T] . 

From now on, we make the additional assumption that our base set of speeds c* satisfies the 
inequalities 

c{<---<c%, (42) 

and we set 

a* = imin{ C * fe+1 -4, 1 < k < N - l} > 0. 



As a consequence of (I25p . (141 p and (|42p . we obtain 

Corollary 4. There exist positive numbers 0:3 < 02 and L3 > £2, depending only on c* , such 
that if ct < a.3 and L > L3, then for any t G [0, T], 

Ofc+l(t) " a k (t) > Ofc+i(0) - o fc (0) + a*i > (L - 1) + (43) 
/or eac/i 1 < A; < N — 1, and 



2-K(i))2>^, (44) 



/or eac/i 1 < k < N. 



E Evolution in time of Q and monotonicity formula 

We define the function Q on [0, T] by 

<?(t) :=G(*7(-, *),«(-,*)), 
where G was introduced in Subsection [Cj We rewrite according to the equality 

JV 

g{t) = E(t) + c \p(t) + £ (4 - 4-i)Qk(t), (45) 

k=2 

where 

tf(t) := £7(T7(.,t),«(.,t)), P(t) := P(r)(;t),v(;t)), 

and 

TV 

Q*(t) := / * fc (a;)p(j7(x,t),u(a;,t))da; = VPj(t). 

j=fe 

In case A = 1 (i.e. there is only one soliton), it follows from the conservation of energy and 
momentum that Q is constant in time. In the general case, we have the following monotonicity 
formula. 

Proposition 5. There exist positive numbers 0:4 < 03 and L4 > L3, depending only on c* , such 
that if a < 04 and L > L^, then 

| (ei(())£o(exp( _M^ 

for any t G [0, T] . 
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As a consequence of the ordering condition (|42|) , we thus obtain 
Corollary 5. If a < 04 and L > L4, then 

'Q{t) )<0 exp' 



cit v w y - v r v i6 

for any t € [0, T] . 

F Uniform control on \\e\\x 

We are now in position to obtain a uniform control on ||e(-,t)||x combining the previous results 
with simple algebraic identities. We divide the analysis in a number of steps, each one providing 
a uniform control on an intermediate quantity. 

Step 1. Control on \Q k (t) - Q fc (0)|. 

Since E and P are conserved by the flow, we obtain from (|45l) . 



Q 3 (0) - Q 3 (t) = 0(0) - Q{t) + — L - 4-i)(Qfc(*) " Qk(0)), (46) 

S' S-i k & 

for any 2 < j < N. On the one hand, we have by (00]) for 0(0), and (J39l) for Q{t), 

^(O)-^(i) < o(|| £ (-,0)||^)+o(|k(-,t)||^)+o(^|c fc (i)-4| 2 )+o( J Lexp(-^-)). (47) 

fc=i 

On the other hand, from Proposition [5] and assumption fj42|) on the speeds, we have 

\ ^(4-4_ 1 )(Q fc (t)-Q fc ( )) <0(exp(-^)). (48) 

Combining (|46p . (j4"T|) and (08]) for the upper bound part, and using once more Proposition [5] for 
the lower bound part, we derive 



N N 



^|Q fc W-g fc (0)|=O(||e(-,0)||i)+O(||e(-,t)||^)+O(^|c fc ( i )-4| 2 )+o(^exp (-^ 

k=2 * k=l 

(49) 

Step 2. Control on \P(Q Ck{t )) ~ P(Qc k (o))\- 
It follows from Proposition [3] that 

P{t) - Q 2 (t) = Pi(t) = P(Q Cl{t) ) + o(\\e(;t)\\x) + °( Lex P(-l^))' ( 5 °) 
and that for 2 < k < N, 

Q k (t) - Q k+1 (t) = P k (t) = P(Q Ck(t) ) + o(\\e(;t)\\x) + O^Lexp ( - ^)). (51) 

Combining (I50p and ()5ip with (j4"9|) and the conservation of the momentum P, we obtain 

N N 

£ - p W Cfc (o))| = o)||i) + e>(lK-,t)||!) + o( E M') " c ^l 2 ) 

k=l k=l (52) 



+ O (l exp 



16 
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Step 3. Control on \ck(t) — c£|. 

Combining (|52p with the fact that 

d 



^- c (P(Qc)) = (2-c 2 )^0, 
and the inequality 



N N N 

Y \c k (t) -4\<J2 h(t) - Cfc (o)| + |c fc (o) -4\, 
it=i fe=i it=i 

we obtain the existence of positive numbers 05 < 04 and L5 > L4, depending only on c*, such 
that if a < 05 and L > L5, then 

AT AT r 



2 |cfc(t) - 4| = O(||e(-,0)|| 2 ) + 0(\\e(;t)\\ 2 ) + 0( £ |c fc (0) - 4|) + OfLexp 



fe=i fc=i 
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N 

Finally, from Proposition [2] we may bound ^ |cfc(0) — c£| by C(||e(-, 0)||x)j and therefore obtain 

A- 1 

Y, h(t) - c%\ = O(|K-,0)||) + 0(\\e(;t)\\ 2 ) +o(Lexp(- ^-)). (53) 



fc=i 

Step 4. Control on \\e(-,t)\\x- 



Combining Corollary [3] for (rj,v) := (i](-,t),v(-,t)), and Corollary [5l we obtain, if a < a.4 and 
L > L 4 , 

A* ^ r 

T \\e(;t)\\ 2 x = o(|k(-,0)|| 2 ,)+o(^|c fc (t)-4| 2 )+o(Lexp(-^-))+o(|K,t)||i). 

fc=i 

It follows therefore from (|53|) that there exist positive numbers ot§ < Q5 and Lq > L5, depending 
only on c* , such that if a < and L > Lq, then 

IK-,*)!! 2 , <o(|K-,0)|| 2 ,)+o(Lexp(-^)). (54) 
G Proof of Theorem [2] completed 

Let (j] ,v ) be as in the statement of Theorem [2j We first impose the condition a* < and 
L* > Lq + 2, so that by continuity of the flow, since L° > L* by assumption, 

T stop :=sup{t>0, s.t. ( V (;s),v(;s)) eU c *(a 6 ,L -2), Vs 6 [0,t]} > 0. 

Notice that T s t op < T max since W c *(ct 6 , L° - 2) C A/"V(M). 
By (|54p we have 

|| e (-,t)|| x = 0(a° + exp(-^)). (55) 



Combining (|55jl with (J53|) and (J54J) yields a positive number .fTg , depending only on c* , such that 
\\( V (-,t),v(;t)) - R c * >a{t) \\ x < K 6 (a° + exp ( - ^)), (56) 
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for any t £ [0,T s t op )- On the other hand, by definition of L°. we have (r](0),v(0)) G U c *(a* , L°). 
In particular, by Proposition [2] 

min (a fc+ i(0) - a fc (0), 1 < fc < JV - 1} > L° - 1, 

so that, by Corollary U we finally have 

mm{a k +i(t)-a k (t), l<k<N-l} = min (a fc+ i(0) - a fc (0), 1 < fe < N - l} > L° - 2, (57) 

for any i G [0,T stop ). 

We therefore additionally impose the condition that a* and L* satisfy 

K e (a* + exp ( - ^p)) < "6, 

so that 

^stop = +00, 

and in particular, ([56|) and (|57|l hold for any t £ [0, +oo). Finally, combining Proposition [4] and 
we obtain 

^K(t)-4|<^ 7 (a° + exp(-^)), 
fc=i 

and the proof of Theorem [2] is completed by setting A* = max{i^6, K-?}. 
1.5 Outline of the paper 

In Section [21 we present the proofs of the results stated in Subsections |Aj [B] and [C] of the 
introduction, which are related to the minimizing properties of solitons and sums of solitons. 
Section [3] is instead devoted to the proofs of the results related to the dynamical properties of 
(jHGPp . and which are stated in Subsections |Pl and [El of the introduction. Finally, for the sake of 
completeness, we provide in Appendix [A] a quantitative version of the implicit function theorem 
which we use in Section [2] in order to define the modulation parameters. 

2 Minimizing properties around a soliton and a sum of solitons 
2.1 Minimizing properties of solitons 

In this subsection we present the proof of Proposition [TJ It is reminiscent from arguments 
developed in [12J. 



Proof of Proposition^ Let us recall that the quadratic form H c is defined in (|20p as 
Equation (Hip writes for r] c as 

8 2 xxVc - (2 - c 2 ) Vc + 3r? c 2 = 0. (2.1) 
Invoking ([5]) and ([2. II) . the expression of H c may be recast as 

H c {e) = L c (e v ) + [ (1 - Vc ) (e v + - e v Y , (2.2) 

JR v ^l 1 — Vc) ' 
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where 



1 f (d x er,r , 1 /" 2-c 2 -6 Vc + 3 V 2 c 2 



The Sturm-Liouville operator £ c associated to the quadratic form L c , namely 



C(s)-=-d( ^ 1 + 2-^-6^ + 3^ 



is a self-adjoint, unbounded operator on L 2 (R), with domain i/ 2 (IR). Since rj c (x) — > as |x| — )• 
+oo by ([5]), it follows from the Weyl criterion that its essential spectrum is equal to [(2 — 
c 2 )/4, +oo). On the other hand, one can translate the invariance with respect to translations of 
equation (jHGPp into the property that the function d x r) c belongs to the kernel of C c . Again in 
view of ([5]), the function d x r] c has exactly one zero. As a consequence, one can infer from standard 
Sturm-Liouville theory (see e.g. [Q\) that the operator C c has exactly one negative eigenvalue 
— Moreover, the associated eigenspace, as well as the kernel of C c , are of dimension one. We 
will denote r]~ an eigenfunction of C c for the eigenvalue — . Notice in particular, that there 
exists a constant < fi + < 2 — c 2 such that 

L c {e v ) >/x + ||^|| 2 2 , (2.3) 

as soon as {e v ,r]~) L 2 = {e v , d x r] c ) L 2 = 0. 

Coming back to the quadratic form H c , we notice that the associated operator 7~L C , namely 

2 

n c (e) := (C c (e v ) + ^^e, + ^4^f», + (1 - Vc)s v ) , 

is a self-adjoint, unbounded operator on L 2 (M) 2 , with domain H 2 (M.) x L 2 (M). Moreover, it 
again follows from the Weyl criterion that its essential spectrum is equal to [(2 — c 2 )/(3 + 
\/T+ 4c 2 ), +oo). In view of © and f|2 . 2[> . we next check that 

Hc { V , ~ 2(l- r]) 2 ) < °' and H c{dxr]c,d x v c ) = 0, 

so that 7l c has at least one negative eigenvalue, and the dimension of its kernel is at least 
one. Assume next that either T~L C owns another negative eigenvalue, or its kernel is at least 
of dimension two. Then, there exists a non-positive direction e = (e v ,e v ) for H c such that 
{£r),V~)L 2 (R) = ^x'?c)l 2 (r) = 0. In view of (|2,2p . this is in contradiction with f|2.3[) . Therefore, 
T-L c has exactly one negative eigenvalue — A~, with eigenfunction \~ := (XrpX^)> while its 
kernel is spanned by d x Q c = (d x r] c ,d x v c ). As a consequence, there exists a positive constant 
A+ < (2 - c 2 )/(3 + y/l + 4c 2 ) such that 

H c (e) > A + ||e||i 2xi2 := A+dl^H^ + || £ „|| 2 2 ), (2.4) 
for any pair e in the closed subspace 

P+ := {(e v ,e v ) G H^R) x L 2 (R), s.t. (s, X ')l^ = (e,d x Q c } L 2 xL i = o}. (2.5) 

The next step in the proof consists in checking that estimate (|2.4p remains available (up to 
a further choice of A + ), when the orthogonality conditions in (|2.5p are replaced by conditions 
([22|) . The proof is reminiscent from We consider the map S(c) := E(Q C ) + cP(Q c ), which 
is well-defined and smooth on (— v2, v2) \ {0} in view of ()18p and (|19p . Using (fTT|) . we compute 

S'(c) = E'(Q c ){d c Q c ) + cP'(Q c )(8 c Q c ) + P(Q C ) = P{Q C ), 
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so that, by (|T8j) . 

S"{c) = P'(Q c ){d c Q c ) = (2-c 2 )i 

On the other hand, taking the derivative with respect to c of the identity E'(Q C ) + cP'(Q c ) = 0, 
we get 

(H c (d c Q c ),e) L 2 xL 2 + P'(Q c )(e) = 0, (2.6) 
for any e 6 A(IR). For e = d c Q c , this gives 

# C (3 C Q C ) = -P'{Q c ){d c Qc) = -(2 - c 2 )5 < 0. (2.7) 

At this stage, we can decompose d c Q c as d c Q c = a\~ + ftd x Q c + r c , with r c G P+. As a 
consequence, 

H c (d c Q c ) = -A_a 2 + H c (r c ), 

so that inequality ([2.7p is equivalent to the fact that a ^ 0, combined with the existence of a 
number < 5 < 1 such that 

tf c (r c ) = ,5A-a 2 . (2.8) 

According to the value of 5, two situations can occur. If 5 = 0, then r c is equal to 0. Since 
a / 0, conditions (|22|) are actually identical to the conditions in (|2.5p . so that f|2.4|) remains 
available under conditions (122 p . If J > 0, we can decompose a pair e which satisfies conditions 
(|22|) . as e := a\~ + ?*, with r £ P c + , and compute 

ff c (e) = -\~a 2 + ^^(r) + ^H c (r). 

In view of (|2.4p . the quadratic form H c is positive on . Therefore, we can apply the Cauchy- 
Schwarz inequality to write 

^ (£) " " A ° + 2tf^) + — W 

At this stage, it follows from (j22|) and (I2.6P that 

{U c {r c ),r) L 2 xL 2 = (n c (d c Q c ),e) L 2 xL 2 + \~aa = \~aa. 
Combining with (|2.8p . we are led to 

H c (e)>^(\-a 2 + H c (r)), 



so that by RIM . 

I — S _ M 2 

flc(e) > — min { A_ ' A+ }!klL 2 xL 2 - 

As a conclusion, up to a possible further choice of the positive constant A + , estimate (|2.4p still 
holds when e satisfies conditions (1221). 



In order to complete the proof of estimate (|2ip . it remains to replace the L 2 -norm of in 
2.411 by an i7 1 -norm. Given a number < r < 1, we write 



_ r I + , _ 1 I (^_ ) + _ 

4ii 1-% V 4 7 R 1 - rj c J 

In view of formulae ([5]) and (j20|) . there exists a positive constant i^T c , possibly depending on c, 
such that 



1 - r? c 



< llrll 2 
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so that, under conditions (12211 . we can take benefit of (j2.4ft to get 

Hc{£) *i I tSt + ((1 - T)A+ - r ^ I'*- 2 - 

4 L Vc 

Since Ty c (0) = (2 — c 2 )/2 = max{|?7 c (x)|, x G M} by ([5]), we deduce that 

Hc{£) - h i {dx£ ^ 2 + ((1 " r)A+ " tK) II^II^xl 2 - 

Choosing r sufficiently small yields estimate (|2ip . Finally, since i? c depends analytically on c 
and the first two eigenvalues of H c are simple, the optimal constant in (|2ip depends continuously 
on c, from which the last statement in Proposition [1] follows. □ 



2.2 Orthogonal decomposition of a chain of solitons 

We now turn to the proof of Proposition [2j One key property in the proof is the exponential 
decay of the functions rj c and v c . As a matter of fact, one can derive from formulae ([5]) that, 
given any integer p, there exists a positive constant K, depending only on p, such that 

£ (\%d*ri c (x)\ + J + * + *\did*v e (x)\) <K(2-c 2 )e W {-(2-c 2 )^\x\). (2.9) 

0<j+k<p 

for any < |c| < \/2 and x G R. A crucial consequence of (|2 . 9[) is the property that two solitons 
with same speed, as well as two sums of solitons with same speeds, can be closed in A(M) only 
if their center(s) of mass are closed. More precisely, we have 

Lemma 2.1. Let c = (c\, . . . , cjy) G Adm^r- Given a positive number 5, there exist two positive 
numbers (3 and M such that, if 

\\Rc,a - -Rc,b||x < &i 

for some positions a £ Posn(M) and b G Posn(M), then, 

N 

^ \a,k - bk\ < 5. 
fe=l 

Proof. The proof is by induction on the integer N. When N = 1, we have 

Step 1. Given a positive number 5, there exists a positive number a such that, if 

\\Qc,a — Qc,b\\x < 

then, 

\a-b\ < 5. 

The proof is by contradiction. Assuming that Step [T] is false, there exist a positive number 5 
and a sequence (a n ) ng N G M N such that 

||Q C (- - O - Qc\\x 0, (2.10) 

as n — > +oo, and simultaneously, \a n \ > 5 for any n G N. If the sequence (tt^)^^^ were 
unbounded, then, up to some subsequence, it would tend to +oo, or — oo. In any case, taking 
the limit n —> +oo into (|2,10p . we would get that 2||Q c ||x = 0, which is not possible. As a 
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consequence, the sequence (a n ) ng N is bounded. Up to a subsequence, it converges to a real 
number a such that \a\ > S, and simultaneously, 

\\Qc(- ~ a) - Q c \\ x = 0. 

This identity implies that Q c (- — na) = Q c for any n G N. Since 0, this ensures, taking the 
limit n — > +oo and invoking (|2.9p . that Q c = 0, which provides the desired contradiction. 

Step 2. End of the proof. 

We now assume that the conclusions of Lemma [2. II are available up to the integer N — 1, and 
that they are not for the integer N. Then, there exist a positive number 5, a sequence (Lj)j e fq 
tending to +oo as j — > +oo, and two sequences (^ n '^)( n ,j)eN 2 an d (^ n '"^)(n,j)eN 2 such that 

a (n,j) G Po SAr (%) and G Posat(%), 

for any n G N, 

||^ Ci n(n,j) - ^c,b("-J) llx(R) °) (2- 11 ) 

as n — )• +oo, and simultaneously, 

AT 

E|a^-6^|>& (2.12) 

k=l 

Without loss of generality, we can moreover assume that a^'^ < b^'^ = 0. 

In this situation, if the sequence (a^'^) n ^fq were unbounded for an integer j, then, up to a 
subsequence, it would tend to — oo as n — > +oo. At this stage, we can invoke (|2 . 9[) to certify the 
existence of a universal constant K such that 

\\R c ,a("' j) ~ RcM n ' j) Wx{M.) >\\R^ a (",j) - R c M n ' j) Wx((-L j /2,+oo)) 

N o J2 ' 



>\\QcJx(8)-K^2^—^( [ exp{-2(2-c 2 k ) l 2\x\)dx) 2 . 

°k \J\x\>LJ2 J 



k=1 - sJ\x\>Lj/2 
for n large enough. In the limit n — > +oo, we deduce combining with (|2.1ip that 

N (2— 2 )j 

WQcJxp) <K^2~ — ex P ( " ( 2 " 4)*Lj), 

k=i Ck 

which is again not possible when j — > +oo. Hence, there exists an integer J such that, given 
any j > J, the sequence (a^^neN is bounded. Up to some subsequence, it converges to a 
non-positive number a J . 

Assume next that a 3 < —Lj/2. For n large enough, we can rely on (12. 9p to write 

l 



N 2 — c 2 i 3L 

+ Ka^> ( X ) " V cM^ ( X ) + u <* (z) I < K ex P ( - ( 2 - c fc) ^ ( x - -^) ) ' 

for any x > —Lj/8. As above, this gives 



\\n ll <ll» /? II i r^f ( 2 - c ^ r J (2-4)^i 

ll^ecjv llx(R) S ||-ft ti0 (n,j) - -K Cj bO,j) ||x(K) + A 2^ eXp 

fc=l Cfc 
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which is not possible at the limit n — > +00, when Lj is large enough. Therefore, we can assume, 
up to a further choice of J, that a? > —Lj/2 for any j > J. 

Invoking once again (|2.9p . we next write 



^2\diVc,a("^( x ) ~ d i%b(^)( x ) ~ d xV CN>a (^){x) + d e x r} CN (x)\ 
£=0 ' N 

N - 1 2-c 2 

+ Ko(».i)( a ~ V cM n ^^ ~ V r»a^^ + V °n( X )\ ^ K Y1 ^ eX P ( - ( 2 - C l) 2 (x- L j)), 

k=l 

for any x > — 3Lj/4, so that 

N—l o, 3 9 1 

lln n II <\\r r II 1 K V ( 2 ~ c fc) 3 n . m f ( 2 ~ c fc)^ L J 

k=l k 

(2.13) 

The expression in the right-hand side of (|2.13p tends to when n — > +00 and Lj — > +00. Up to 
a further choice of J, one can invoke Step [T] to show that 



,0»| 



< 4, (2-14) 



for any j > J and any n > rij. Moreover, one can rephrase (|2.13p as 
N-l 

(Q (n.i) — Q , fn.il ) 

X(R) 



^ ^ Q c k ,a^ ) ~ Q c k A n ' j) ) 
k=l 

" (2-c?)t / (2-cL)iL 
< 2\\R tia{n ,i) - R cMnJ) \\ xm + KY,- — exp ( - ^ f— 



(2.15) 



k=l 



Once again, the expression in the right-hand side of (|2.15p tends to when n —> +00 and 
Lj — > +00. Since Lemma [2~T1 is true for the integer N — 1, we deduce that 

N-l x 

£| a y.»)_ 6? ,»),^ 



k=l 



S 4' 



for any j and n large enough. Combining with (|2.14p . we obtain a contradiction with (|2.12|) . 
Lemma [2.11 follows by induction on the integer N . □ 

Estimate (|2.9p and Lemma 12.11 are enough to consider the proof of Proposition [2J 

Proof of Proposition^ The proof is reminiscent from [16j. It relies on the quantified version of 
the implicit function theorem provided by Proposition IA.1I We consider the map S defined by 

E(rj,v,a,b) := ({e,d x Q aubl ) L 2, . . . , {e 

b N ) tfiP' {QaiM)( £ )i ■ ■ ■ ,P'(Qa N ,b N )( £ ))- 

In order to simplify the notation, we have set here e := (77, v) — R a ,b- The map H is well-defined 
and smooth from x Adm^r x 1H N to M. 2N . Moreover, it fulfills the assumptions of Proposition 

I A. II so that we can state 
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Step 1. Let < r < 1 and set 

Adm7v(r) := {c G (-V2, V2) N , s.t. /i c > r and z^ c > r}. 

There exist positive numbers 5, A and M, depending only on t, such that, given any (c, a) G 
Adni7v(T) x Posat(M) ; there exists a map 7 c n G C^B^a,*), AdniAr(r/2) x 1^) suc/i i/iaf, 
given any (rj,v) G B(R Cta ,5), (a, b) = 7 c , a ( ? ? ) ' u ) * s ^ e unique solution in B((c,a),A5) of the 
equation 

S(ri,v,(T, b) = 0. 

Moreover, the map 7 Ci0 is Lipschitz on B(R Cta ,5), with Lipschitz constant at most A. 
For c G Adm at(t) and a G M. N , we check that 

E(R ca ,c,a) = 0. 



Similarly, we compute 



do-jE.f.(Rc,a, c, fl) — —(d c Q C:jtaj ,d x Q Cktak )i2, 
dbjEkiRca, c, a) = (d x Q c - a -,d x Q Cktak ) L 2, 



as well as 



^6,-3iV+fc(-Rc,a! c ) a ) = P'(Qc k ,a k )(d x Q c . aj ), 



for any 1 < j,k < N. When j = fc, we infer from formulae ([5]) that 

d ak E k (R C:a , c, o) = d ak E N+k (R CA , c, a) = 0, (2.16) 



and 



whereas by ([18 



afc S fc (ii Cia ,c,a) = HSxQcJIia = (2 Q Cfc)2 > 



a^S^^.ca) = -l(p(Q c )) |c=Cfe = -(2-4)1 < 0. 

In particular, the diagonal matrix A c with the same diagonal elements as d as bE(R t a , c, o) is a 
continuous isomorphism from M. 2N to M? N , with operator norm bounded from below by r 3 /3. 



On the other hand, when j ^ k, it follows from (|2,9p that 

| (9 c Q C j )Q . , dxQc^^a^ ) L 2 | "I - | {&xQcj,aj j dxQck^a^ )l 2 | 

+ |P , (Q Cft , a J(a c Q Cj , aj )| + |P / (Q Cfc ,a fc )(S a! Q« y ,a J )| 
1 

1 

where we have set i/^jj. := min{(2 — c^) 2 ,(2 — c 2 ) 2 }. Assuming that a G Posjv(L) for some 
positive number L, there exists a positive constant K T , depending only on r, such that 

I^.Hfe^c^c,^! + \d aj E k (R c>a ,c,a)\ 

+ \d aj E N+k (Rc, a , c, a)\ + \d aj E N+k (R^ a ,c,a)\ < K T (L + l)exp ( - u t L). 



< K 2 (2-c|) 2 (2-4) 2 (l + ^- T ) / exp(-i/ ijfc (|a!-o^| + |x-a fc |))da; 
= K 2 (2 - c])5(2 - cl)Ul + -^3) (\aj - a k \ + —) exp ( - u jtJt \aj - a k \), 
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Combining with (|2,16p . we deduce when L is large enough, that 

d a!b E(R^ a , c, o) := A c (id + H(c, a)) , 

where the operator norm of H(c, a) is less than 1/2. 
We now turn to the differential 

d v ,vE(Rc,a, c,o)(/,ff) := (i(f,g),d x Q cltai ) L 2,. . . ,((f,g),d x Q CN 

An) L 2 i 

P'(Q cuai )(f,g),...,P'(Qc N ,aM,g))- 

Since the operator norm of A~ l is bounded by 3/r 3 , we infer from (|2.9p that cL w 5(i? c , c, a) 
may be written as 

where T C (J is a continuous linear mapping from X(R) to with operator norm depending only 
on r. 

Finally, again by (|2.9|) . the operator norm of the second order differential d 2 H(r/, v, a, b) is 
bounded by a constant K T , depending only on r, when (rj, v, a, b) G X(R) x Adm7v(r/2) x R N . It 
then remains to notice that Assumption (iv) of Proposition [Aj] is satisfied when U = Admjv(T) 
and V = Admjv(T/2), and to apply Proposition IA.1I to the map S, in order to obtain the 
statements in Step[TJ 

Step 2. End of the proof. 

Let r = min{/i*/2, u*/2}. We denote 8\, Ai and Mi, the constants provided by Step [TJ 
and /?i, the number provided by Lemma [2.11 for 8 = Ai<5i/3. We set a\ := min{<5i/3, /3i/4} and 
Li := Ml. When (r],v) G U c *(a\, L\), there exists b G Postv(-^i) such that (rj,v) G B(R ( -*^, 2a%). 
Since 2«i < #i and Li = Mi, we deduce from Step [1] that the numbers c and a given by 

(c,a) = j c * ib (r),v), 

are well-defined, so that we can set <t(r], v) = c and 2l(r/, u) = o. We claim that c and a do not 
depend on the choice of b G Posat(Li) such that (77, v) G B(R C *^, 2a\). As a consequence, the 
maps (£ and 21 are well-defined on U c *(ai,Li) with values in AdmAr(r/2) C Admjv, respectively, 
R N . 

Indeed, given another choice of b<i G Posjv(Li) such that (rj, v) G B(R C * & 2 , 2ai), we have 

\\Rc*,b ~ Rc*m\\x < 4ai < /Si, 

so that, by Lemma |2~T1 

|b-M<^- (2-17) 

On the other hand, the map "f c *,b 1S Lipschitz on B(R c *j,, 2a%), with Lipschitz constant at most 
Ai. Hence, 

|(c,c) - (c*,b)| = |7e,b(»7,u) - 7c-,b(-Rc-,b)| < Aiai < — 
Combining with (|2,17p . we obtain 

|(c,o) - (c*,b 2 )| < Ai8i. 

Since E(r],v, c, a) = 0, we deduce from Step Q] that (c, a) = 7c*,b 2 ( r 7> u )) so that c and do not 
depend on the choice of b G Postv(-^i) such that (rj, v) G i?(i? c * 5 [,, 2a\). 
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Concerning the smoothness of £ and 21, we consider (n,w) S W c *(ai,Li) such that 

HO, w) - (v,v)\\x < ai. 
Given b E Posat(Li) such that (77,71) S B(R c * t b,ai), we deduce that 

\\(n,w) - R c *,b\\x < 2«i, 

so that, by definition, 

(<£(r],v),$l(r),v)) = -y t *,b(v,v) and (t(n,w), 2t(/i, w)) = 7 c * ) &(/x,w). 

At this stage, recall that the map is of class C 1 on B(R C * &, 2ai). As a consequence, the 
functions (£ and 21 are also of class C 1 on B((rj, v), ai) nU c *(az\,Li), therefore, on 1A t *{oti,L\). 

We next consider the proof of (|25|) . When (rj,v) £ -B(i? c * jQ * , a), with 0* S Posat(L), a < ai 
and L > L\, the values of C(rj,v) and 21(77,7;) are equal to 

(£(t/,v),21(»7,v)) = 7c*, a * 
In view of the Lipschitz continuity on B(R (: *^ a *, a) of the map 7 C *, *) we infer that 

|c-c*| + |o-a*| < AiUv,v) -Rf,a'\\x < Aia. (2.18) 
On the other hand, letting s = (77, t> ) — -R c ,a> we have 

Ikllx < \\(rj,v) -R c *, a *\\x + \\Rc*, a * -Rc, a \\x < a+ \\R c *, a * - Rc,a\\x- (2-19) 
In view of formulae ([5]), there exists a universal constant K such that 

\\R c *, a * - R c ,a\\x < / 3 CT iVK(c*-c),a+t(a*-a)( c * - c) + <9*i? c _K(c*-c),a+t(a*-a) ( a * - 0) 
J 

<tf|(c*,a*)-(c,a)|. 
Combining with (|2.18p and f|2. 19[) . we are led to 

\\e\\x < (l + ATAiV 



dt 

x 



In view of (|2.18p . it is sufficient to set K\ := 1 + {K + l)Ai in order to derive ()25|) . 

Finally, conditions (|24p are direct consequences of the definitions of the maps 7 C * &. This 
completes the proof of Proposition [2j □ 

Corollary [1] is a direct consequence of Proposition [2l 

Proof of CorollaryUl Let (77,1") G U c *(a\,L), with L > Li, and set (c, a) := (£(77, v), 2l(?7, w)) as 
in the proof of Proposition [21 In view of (|25|) . one can decrease a\ so that 

1 *i ^ • r/v 1 

c — c < mm <^ , — — >, 

' ' " I 2 8V2) 

which is enough to obtain (I28p and (|29D . 

One can next increase L\ so that L\ > 2, and decrease again a\ so that K\oi\ < 1/2. In this 
case, we can deduce from j25j) that 

a — a < -, 
1 2 
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for any a* G Posn(L) such that 

\\{r],v) - R c *,a*\\x < a. 
It follows that a belongs to Posn(L — 1) as mentioned in (\27\) . 
Finally, one can infer from formula ([5]) that 



mm 



[l-Vc(x), xGK} = l-7/ c (0) = ^ 



for any c G (— y2, \/2). Combining with the exponential decay of the functions r/ c provided by 
(|2.9I) . we deduce that we can enlarge again L\ so that 

minjl -7/ c ,a(x), x G m} > (2.20) 

when c satisfies (|29p and o G Posat(Li — 1). On the other hand, we can a last time decrease a\ 
so that, combining (|25p with the Sobolev embedding theorem, we have 

II* II < & 

H e 77l|Z,°°(R) S -g~- 

Since r] = 7y Cj9 + e„, this is sufficient to guarantee that r] satisfies (|26p . so that the pair (77, v) 
belongs to AAV(K). This completes the proof of Corollary [TJ □ 

2.3 Almost minimizing properties close to a sum of solitons 

The main goal of this subsection is to establish the almost minimizing properties close to a sum 
of solitons stated in Subsection 0°f the introduction. When establishing the estimates of E and 
Pk given by (|33[) and (|34p . we will make use of the following elementary inequality to bound 
various interaction terms using (|2.9p . 

Lemma 2.2. Let (a, b) G 1R 2 ; with a < b, (u a , v b ) G (0, +oo) 2 7 and set y^ := max{±y, 0}. Then, 

2 ~ 

exp (-v a (--a) + ) exp (-u b (--by) < ( — : T + b- a ) P exp ( - mm{u a , v b }(b- a)) , 

for any 1 < p < +00. 



We will also use the following pointwise estimates on the functions $kk+l an d ^k- 
Lemma 2.3. Let 1 < k < N and i£R. We have 

Li- 



®k(x) < exp ^ - 2t[\x - a k \ 
and 

|1 - *fe(a?)| < exp ( - 2r(x - a fc + -^j ) + exp ( - 2t(x - a fc - -j- 
Similarly, for 1 < k < N — 1, 

Li 



$q,i(x) < exp ^ - 2t[x - ai + ^ 
$fc,fc+i(a;) < exp ( - 2t(x - a*. - ^) ) exp ( - 2r(x - a k+1 + ^ 
$Ar,Ar+i(s) < exp ( - 2r^x - a N - — ^ 
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while, for 2 < k < N - 1, 

- *2(jc) < exp - — (x 

I'fcO) - *fc+i(a;) < exp ( - ^- (a; - ak ~ l + ak ^j ^ e xp ( - ^ (x - afc ~^ + afc 
*iv(x) - Wjv+i(x) < exp - — (^x ; 

and 

|1 - *i(x) + * 2 (x)| < exp (- ^ (x - 

|1 - *,(*) + < exp ( - ^(x - + exp ( - - 

|1 _ + ^ +1 ( X )| < exp ( - V f (x - a »-i+ a » 

Proof. The estimates in Lemma [2.31 follow from the inequality 

|l — sign(x)(x)| < 2exp(— 2|x|), 
which holds for any i£l, □ 

We are now in position to present the 

Proof of Proposition [5j We begin with (|33[) . We expand the energy E according to the Taylor 
rule by 

E{ri, v) = E{R c>a + e) = E(R CA ) + E'(R c>a )(e) + ^"(i? c , a )( £ , e) + Tl( V , v), 



where 



(da^ca) 2 I,, x 2 , 1 2 



E (Rc, a )(e) = / — {- + — ( + (1 - J7 Cl a)ucae« + ^(»7c,o - «c,a) e *j , (2.22) 

£ (Rc,a){e, e) - ^ ^^(^^ + 2(1 _ ^ + 1(73^)3 " + (l " Vc,a)s v + y j , 

(2.23) 

and 

^ , _ 1 f ( {dxVc, a ) 2 £* | 2fl E 7; I , a Q 1! e t; e;) (9^)% , 2 \ 

8 7 M l(l-r ? )(l-r ?Cj0 )3 + (l-r / )(l- ?7C)0 )2 + (l- ?7 )(l-r 7cja ) 

Following Corollary [H we may assume that (j26j) . (|27|) . (j28[) and (l29|) . as well as estimate (|2.20p . 
hold. From (j2~9|) and ([25)1. we obtain 

\ric k , ak (x)\ + IflUifc.^Cx)! + \v Ck , ak (x)\ = O ( exp ( - ZV|X ~ afcl )), (2.24) 

for any x S M and any 1 < k < N. From ((2J), ff2~20]) . (l2~24")l . and the expression (|2T2Tj) for 
a ), we compute using Lemma [2^2] with j> — 1. 

iV 

#(i?c,a) = E ) + ° ( L eX P ( - ^) ) • ( 2 ' 25 ) 



k=l 
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Similarly, expression (|2,22p for E 1 may be simplified using the Holder inequality and Lemma |2~2 
with p = 2 as 



N T 



E\R tA ){e) = Y,E'(Q Ck , ak )(e) + oilA exp 



\ e \\x 

\ \ L / 

k=l 

Since \\e\\x < K±ai by (|25p and since 

E'(Q Ck , ak )(e) = E'{Q Ck . ak ){e) + c k P'(Q Ck , ak )(e) = 0, 
by (HB and ([Ml), we deduce that 

E'{R CA ){e) = 0(h\ exp ( - ^)). (2.26) 

From (|26p . (|2.20p and (I2.24p . as well as the Sobolev embedding theorem, we also infer that 

K( V ,v) = 0{\\e\\ x ). (2.27) 

We finally turn to the term E"(R Cta )(e,e) which requires somewhat more care. We rewrite 
expression (|2.23p as 

E"(R ca ){e,e) := [ J(R c , a ,e), 

JR 

and we decompose according to the partition of unity (|30p as 



„ N N 

/ J{R C:Q ,e) = Y^ / J(Rc,a,e)^k + Yl / J{Rc,a,e)^ k ,k+i- 

JR , — -, JR JR 



k=l ^ k=0 



Invoking once more the decay estimates (|2.24p . we obtain from Lemma 12.21 with p = +oo, and 
Lemma 12.31 

TLl \\\ 2 x). (2-28) 



/ J(R C}Sl ,£)$ k = / J(Q Cfc , afc ,e)$ fc + o(exp 

JR JR v 

Taking into account the formula 



d x (<!>le r] )=<S>ld x e v + ^e v , 
we obtain, after a change of variable corresponding to a translation by a k , 



/ J(Q Ck ,a k ,e)$k 

JR 



32<D fc (l-77 Cfc , a J 8(1-7?^) 8(1 - r?, 



c fci a fc J 



(2.29) 



Since \d x <S> k \ = 0(r) and \d x $ k \/\<b k \ = O(r), we deduce from (l2~2"9l> . 

/ J(Q Ck ,a k ,e)^ k = E , \Q Ck )(e k ,e k ) + 0(r\\e\\ x ). (2.30) 

JR 

In a similar manner, we write using Lemma 12.21 with p = +oo, and Lemma [2.31 
/ J(i? c , a ,£)$ fc>fc+ i = / J(0,e)$ k . k+l + C(exp 

JR JR v 
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IX J. 



and then 



/ J(0,e)<S> kik+1 = [ < 7(O J e fc , JW .i) + 0(r||e|ft) = E"(0)(e ktk+1 ,e kik+1 ) + 0(r\\ef x ) 
Jr Jr 

Estimate ([33]) follows combining with ([2~25t ([2~26t (l2~2Tl) . (pT28|) and ([2~30l) . 

We turn now to (|34p . We expand the functional P k , which is exactly quadratic, as 



P k (v,v) = P k (R c , a + e) = P k (R c , a ) + P' k (R c , a )(e) + ~P£(i? C)0 )( £ , s), 



where 



Pk(Rc,a) = l [ ^c,a(*fc-*fc+l), (2-31) 

1 Jr 

P k (R c ,a)(e) = \ [ric, a e v + v c>a e v ) (* fc - , (2.32) 

z JR 

i^'(i2 c , tt )(e ) e)= / e^e„(* fc - * fe+ i). (2.33) 
Jr 

Invoking here again ([2.24|) . Lemma 12.21 with p = 1, and Lemma |2.3| we estimate (I2.3ip as 

N j 
Pk(Rc,a) = Yl PkiQc^) + O^Lexp ( - 
i=i 

For j 7^ A;, we obtain from Lemma [2.21 with p = 1, and Lemma |2.3[ 



P fc (Q c „a i ) = 0^exp v |( . 

while 

i\(Oc h ,«J = P(Qc k ) + ©(Lexp ( - ^ 
Concerning (12.321) . we obtain similarly from Lemma [2.21 with p = 2, and Lemma [273 

^(iU(e) = m*,«J0O + 0{Lh exp ( - ^)) = C?(l§ exp ( - ^ 

where we have used (|24l) for the last equality. 
Finally, we decompose the integral in (|2.33|) as 

„ N r N „ 

/ = V / -^ fc+ i)$j + V / £^(# fc -* fc+ i)$ i)i+ i. 

JR j =1 JR J=0 JR 

When j ^ k, we obtain 

J e v £ v (y k - V k +i)®j = £>(exp (- ^ i )||e||x), 

while for j = k, 

/ e^e„(*fc - = / (e^$|)(£ v $|) + / e v e v <Z> k (^ k - ^ k+1 - l) 

JR JR 

= P"(Qc k )(e k ,e k ) + Of exp f - T%L) \\e\\ 2 x 
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Also, when j ^ k and j ^ k — 1, we obtain 

£ v e v ($ k - * fc+ i)* iij+ i = ©(exp ( - ^) 
The remaining two terms (corresponding to j = k and j = A; — 1) are simply grouped together as 

/ S^v^k — ^k+l) {^k-l.k + ®k,k+l) = P'k (®)( £ k-l,k, £ k-l,k) + ffe (0)( e fc,fc+l> £ft,fc+l)' 
./R 

The proof of estimate (I34p . and therefore of Proposition [3J is completed. □ 

As a consequence of Proposition we derive 
Proof of Corollary \^ We write 

E(Q Ck ) = E{Q cl ) + E'(Q c *)(Q Ck - Q c% ) + O(\o k - c%\ 2 ) , 

and 

P{Q Ck ) = P(Q ci ) + P'(Qct)(Q Ck - Qc* h ) + O{\o k - c* k \ 2 ) . 

In view of Proposition [3j the conclusion follows from the identity 

E'(Q cl ) + 4P'(Q cl )=0, 
and from the inequalities (in order to recover the operators H Ck and Hq ) 

\(4-c k )P"(Q Ck )(e k ,e k )\ < \4-c k \ [ \s v \\e v \ < O{\o k - c%\ 2 ) + 0(\\e\\ x ), 

Jr 



and 



|(4 - c k )P^0)(e j>j+1 ,s j)j+1 )\ < |4 - c k \ [ \e v \\e v \ < 0{\c k - 4| 2 ) + 0{\\e\\ x ) 

Jr 



□ 



We next prove Lemma [TJ 



Proof of Lemma[Jl Let 1 < k < N. Under the assumptions of Corollary [2l the pair e satisfies 
the orthogonality conditions (|24p . In particular, we can write 

{£k,d x Q Ck ) L 2 = / (s,d x Q Ck , ak ) R 2<S> k = / (e,d x Q Ch)ak )&i($k-l), 
Jr Jr 

and 

P'{Qc k )(e k ) = - / {£,(v Ck , ak ,ri Ck>ak )) R 2<£> k = - / (e,(v Ck>ak ,ri Ckyak )) R 2($ k - l). 



Arguing as in the proof of Proposition [3l we infer that 

J \e\{\d x Q Ck ,a k \ + \(vc k ,a k ,Vc k ,a k )\)\$k - l| = O (^f exp ( - ^) ||e|| L 2 

so that ^ 

(£k,d x Q Ck ) L 2 = 0[L{ exp ( ^-jllellis), 
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and 



P'(Q Ck )(e k ) = 0\L\ exp 



\£\\l 2 • 



We next write the following orthogonal decomposition of e k , 

(e k ,d x Q Ck ) L 2 (e k ,(v Ck ,r] Ck )) L 2 



£k ■-- 



\\d x Q 



c k \\L 2 



-u k 



{v Ck ,ri Ck )\\ L 2 



Vk + r k , 



(2.34) 



(2.35) 



in which u k := d x Q Ck /\\d x Q Ck \\ L 2 , w fc := (w Cfc , 77 Cfe )/|| (^c fe , *?c fe )[|i2, and r k satisfies the orthogonal- 
ity conditions (|22|) . In view of Proposition [TJ we have 

ff Cfc (r fe )>A c Jr fc ||^. 



On the other hand, since d x Q Ck and hence u k belongs to the kernel of the quadratic form H t 



<-k ■ 



we can write 



n , , ;/ , , , „ P'(Qc k )(e k ) m , , . , {P'{Q Ck ){e k )) u , , 
H Ck {e k ) = H Ck {r k ) + 2— — — {n Ck (v k ),e k ) L 2 + — — H Ck (v k ). 



ivc k ,r) Ck )\\ L 2 



L 2 



It follows from (|2~33|) that 



H Ck (e k )>A Ck \\r k \\ 2 x + Lf exp 



L 2 



(2.36) 



Next, we similarly deduce from (|2.35p that 



\r k \\x ~ IkfclUl < 



{^k,d x Q Ck ) L 2 



dxQc k \\l 2 
and therefore from (|2.36p . 



P'(QcJ(^) 



x = C£fexp 



l e IU 2 



H Ck (e k ) >^\\ek\\ 2 x + 0(l\ exp ( - T J± 



\4\ 



(2.37) 



We now turn to the terms HQ{e k)k+ \). For 1 < A; < TV — 1, we write 

[c fe P£(0) + c fc+1 P^ + i(°)] (E*,fc+i,efc,*+i) 

= / c k (^k - + Cfc+i^fe+i - *fc +2 ) (efe,fc+i) (efe,fc+i) ■ 



Since 
we obtain 



|c fc (*fc - + c fc+ i(^fc+i - ^fc+2)| < max{|c fc |, |c fc+ i|}, 



H$(e k>k+1 ) > [E"(Q) - max{|c fc |, |cfc+i|}P"(0)] (£ k)k+ i,e k)k+1 ). 
On the other hand, the quadratic form E"(0) — v2P"(0) is non- negative, so that 
max{|c fc |,|c fc+ i|}\ „ max{|c fc |,|c fe+ i|} 



Ho(e k , k +i) > ( 1 
Similarly, 



V2 



H° (s 0>1 ) > (l - ^)P"(0)Ki,eo,i) >\(l 



\ci\ 
V2 



V2 



lko,i||x, 



lkfe,fc+illx- 

(2.38) 
(2.39) 
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and 



H»(e N>N+1 ) > (l - M)E"(0)( £W1)£W1 ) > \{l - M)\\ £NtN+l f x . (2.40) 



In view of (l2~3T)l . (l2~38l) . ([2~39l) and (f2~30|) . we set 

' A c fc 1/, |Cfc| 



A * = mm ii"'4l 1 -^J' 1 ^^ iV i' 

so that (|35p and (|36[) follow. Notice that Corollary [1] and the last statement of Proposition [1] 
ensure that A* only depends on c*. 

In order to complete the proof of Lemma Q] we now derive ([37p . In view of definitions (|3 1 1) 
and (|32p . we may write 

N N 1 

EiNi* = E/ ((^*i)) 2 + M4+ e *))> 
/c=i fc=i' /R 

and 

N N 1 

Ell e M+illx = £ / ((Me^l k+ i)) 2 + ^k,k + i{e 2 v + el)). 

k=0 k=0 jR 

l l 

Expanding of the derivatives d x (e rj ^^.) and ^(e^^ and summing, we are led to 

N N N N 



£iniI+£ik*+iIIx = / (E $fc + E^ fc + i )( (ax£r ' )2 +£ v +s v) 

k=l k=0 R k=l k=0 

r N N N r fa if> \z N r 

+ /„ * (I> + E *^)^ + E / R ^4 + E /, 



fc=l fc=0 fc=l"^ * fc=0' 



Since 

N N 

and <I>fc > on 1R (so that the last two sums in the previous identity are non- negative) , we finally 
deduce that 

N N . 

£Niix+£N,*+ii&> / (( d * £ v) 2 + 4 + £ v) = \\ £ \\x- 

k=l k=0 Jr 

This completes the proof of Lemma [TJ □ 

Proof of Corollary [5j Inequality (|39p follows from Corollary [21 Lemma [1] and inequality (|38p . 
while inequality (|40p follows from Corollary [2] and the fact that H Ck and Hq are continuous 
quadratic forms, with continuous bounds depending only on c*, in view of Corollary [U □ 



3 Dynamical properties 

In this section, we focus on the evolution in time under ([HGPp of various quantities introduced 
so far. In particular, we present the proofs of Propositions [4] and [5l 
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3.1 Dynamics of the modulation parameters 

Proof of Proposition^ We first derive the conclusions of Proposition U when the pair (r/,v) is 
a smooth solution to ()HGP[) . for which we are allowed to compute the differential equations 
satisfied by the modulation parameters a k (t) and c k (t). We next complete the proof using a 
density argument. 

Step 1. Regularity and estimates for smooth solutions. 

In this step, we make the further assumption that (rj,v) G C°([0, T], X 3 (M)), from which 
it follows by equation (IHGPp that (r),v) G C 1 ([0, T], X(M)), and then from Proposition [2] that 
c € ^([O^j.Admjv) and a G ^([O.T],^). We derive from dHGPjl that e := (er,,e v ) is solution 
to the equations 



d t e v =2d x (e v - e v v iA - e v r} iA - e n e v ) - y~) d x (r) Cktak v Ci >a . ) 

fc=i jy/c 

JV JV 



(3.1) 



fe=l fc=l 

and 

JV JV 



fe=l j^k k=l 

JV 



(l-Vc k ,a k )' 

a2 ( d x r/ t<a + d x e v \,\^ a f 3(d x n Ck , ak ) 2 \ a ( 3(d x ri c , a + d x e v ) 2 \ , . 
X \2(l-V, a -e v )) + h W~ri Ck ,a k ?) ~ 9x ^ 4(1 - n c , a - ) ^ 

JV JV 

C 'kdcV Ck ,a k + J] (4 - C k )d x V Ck>ak . 



k=l k=l 



Here and in the sequel, we have dropped the explicit mention to the dependence in t when this 
does not lead to a confusion. 



We next differentiate the orthogonality conditions (j22|) to obtain 

j t {e{-,t),d x Q Ck{t)) a k{t) ) L 2 = ^P'(Q Ck (t),a k (t))(£(-,t)) = 0, 

for any 1 < k < N. Expanding the previous identities according to the chain rule and combining 
the resulting equations with (|3.ip and fj3.2j) yields the system 



M « („<( t ) - c( t )) = »(')• < 3 - 3 > 
where the coefficients (M k j(t))i<j tk <2N of the 2iV-dimensional matrix M(t) are given by 



M hj (t) := P'(Q Ck ,a k )(d c Q Cj , aj ) - 6 jik P'(d c Q Ck>ak )(e), 
M kij+N (t) := -P'(Q Ck , ak )(d x Q CjA] ) + 6 j>k P'(d x Q Ck , ak )(e), 
M k+ Nj(t) := (d c Qc J ,a j ,d x Q Ck: a k )L 2 — 3j,k(£,d c d x Q Ck , ak ) L 2, 
M k+ N,j+N(t) ■= -(d x Q Cjtaj ,d x Q Ck , ak ) L 2 + 5 jtk (£,dlQ Ckt a k ) L 2, 
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and the components ($ k (t))i<k<2N of the 2N- vector $(t) are given by 

®k(t) := - 2(e v - E v v c>a - £ v rj Cj a - e v e v , d x v Ck)Clk ) L 2 - (e v - 2v^ a e v - (e v ) 2 , d x rj CkiClk ) L i 

- c k P (d x Q Ck , ah )(e) - { 2{l _^_ £ii) ~ 1, 2(1-^.)' 

/ 3(a a T/ C|0 + a a! g t> ) 2 ^ sjpxVcj^ R n \ (3.4) 



\ 4(1 - ^ c , a - e„) 2 ^ 4(1 - r, Cj>aj ? ' 

N N 



j=l 1=1 



+ 5^5^ i{Vcj,aj V ci,ap d x V Ckt a k ) ip. + {v c ^ a .V Cuav d x r] Ck ^ k ) L 2 J , 



$ k+N (t) :=2{e n v Cfi + e v 7] Cta + e v e v - e v , d%q Ck>a , k ) L 2 - {e^ - 2v^ a e v - (e v ) 2 , d x v Ck>ak ) L 2 



C k {£,0 x y Ck , ak)L 2 ^ x 2^2(1 



■j i a j o3 

5 °x v c k ,a k 



2(1 - r? Cia - e„) ^ 2(1 - ,,,.,.„, )' * "«'-/ ^ 
/ Kdx-n^ + dxSrj) 2 3(^7 ?Cjiaj .) 2 \ (3.5) 



/ <i{O x r] C:a + O x £rir \ - ■!{OxVc J ,a j ) 2 \ 

\ 4(1 - ^ - e,) 2 A. 4(1 - 7?c ,, ai ) 2 ' ^ C&,afc / 



L 2 



AT AT 



,a,j v ci,aii ,a k )l 2 



3=1 1=1 

We next write the matrix M(t) under the form Af(t) := -D(i) + H(t), where -D(t) refers to the 
diagonal matrix with diagonal coefficients 

D k>k (t) := P'(Q Ck )(d c Q Ck ) = (2 - c|)i / 0, 

(2- c 2 )i 

Dk+N,k+N(t) := ||<9;cQcJlz,2 = / 0, 

where 1 < < iV. In particular, -D(i) is invertible, and its inverse is uniformly bounded with 
respect to t in view of (I28|) . In order to estimate the matrix H(t), we first invoke the identities 

P'(Qc k ,a k )(d x Qc k ,a k ) = (dcQc k ,a k ,d x Q Ck: a k ) L 2 = 0- 

Since the remaining terms involve either e or distinct solitons, we can rely on computations 
similar to the ones in the proof of Proposition [3] to get 



\H kjj (t)\ =0(||e|| L2 )+C>(Le X p( 



2 

Hence, there exists 0.2 < ol\ and L2 > L\ such that for a < a.2 and L > L2, the operator norm 
of the matrix D(t) -1 H(t) is less than 1/2, so that M(t) is invertible. Moreover, the operator 
norms of the matrices M{t)~ l are uniformly bounded with respect to t. 

At this stage, we have proved that 

N 

H(t) - cm + \c' k (t)\) = o(n*(*)[| Haw ), 

k=l 

for any t £ [0,T]. Inspection of $(t) shows that, similar to H k j above, it only contains terms 
involving either e or distinct solitons. Therefore we obtain 

N 

E (K(*) - c*(*)| + I4(*)l) = 0(||e|| La ) + o(Lexp ^ ' V " 
fc=i 
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which corresponds to (|4ip . 
Step 2. Density argument. 

We consider a sequence (t]q,Vq) in AAV 3 (M) such that 

as n — > +00, and we denote (n n ,v n ) the corresponding solutions to (jHGPp . It follows from 
Theorem [1] that the solutions (rj n ,v n ) are well-defined on [0, T] for n sufficiently large and that 

( v n ,v n )^( V ,v) in C°([0,T],X(R)), (3.6) 

as n — > +00. In particular, for n large enough, the compactness of the segment [0, T] ensures 
that (n n (■ , t) , v n (■ , t)) belongs to £Y c *(a,L) for any t £ [0, T]. As a consequence, one can apply 
Step [1] to (n n ,v n ). This yields maps c n and a n of class C 1 on [0,T] such that the orthogonality 
conditions (|2'2p are satisfied by the pair e n := (e",e") = (n n ,v n ) — R c n a n. Moreover, combining 
the Lipschitz continuity of the maps £ and 21 on some tubular neighborhood of the compact set 
{(ri(t),v(t)), t G [0,7]} with QSty . we have 

(c"(t),a»(t))->(c(i),a(t)) in C°([0,T],IR 2Ar ), (3.7) 

as n — )• +00. Recall that, by formulae ([5]), the maps (c, a) — > d"dxQ c (' — a ) are continuous from 
((-a/2, \/2) \{0}) x R into L 2 (M) 2 for any (a, /3) E N 2 . Combined with f[577]) . this shows that the 
matrices M n {t) converge towards the matrices M(t), uniformly with respect to t, as n — > +00. 
Since the map A 1— > A^ 1 is continuous from QCn(R) to 7Wat(M), the inverses (M n (i)) -1 also 
converge towards the inverses M(t)~ l uniformly with respect to t G [0, T]. 

We next turn to the vectors $ n (f). In view of (|3,6p and ()3.7p . we can write 

e n = (r ) n ,v n )-R c n ia n^( v ,v)-R Cta = e in C°([0, T],X(M)), (3.8) 

as n — > +00. On the other hand, in view of formulae ([5|), the derivatives d^dxQc have a 
vanishing limit at infinity for any (a, /?) G N 2 . Therefore, the maps — a) are 

not only continuous for the L 2 -norm, but also for the uniform one. Applying (|3.7p and (|3.8p to 
expressions ()3.4p and (|3.5p . and invoking the Sobolev embedding theorem of i? 1 (M) into C (M) 
when necessary, it follows that 

in C°([0,T],M. 2N ), 
as n — > +00. Coming back to (|3.3p . we are finally led to 

U'w-%) M(t)-1 * (t) - 

as n — >■ +00, the convergence being uniform with respect to t 6 [0, T]. In view of (|3.7p . this 
shows that the maps c and are actually of class C 1 on [0, T]. Moreover, they are given by 

U) - c( t )) - M(t)_1 * (t) ' 
for any t G [0, T], so that (|41[) follows as in Step[TJ □ 

We are now in position to provide the 
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Proof of Corollary^ When (rj(-,t),v(-,t)) G U c *(a,L) for any t G [0, T], we can combine esti- 
mates (1251) and (l4"T1) to obtain 



N N N T 

E Wk® - 4| < E K*(*) " c ^)l + E M') " 4| = <%) + ©((iexp ( - ^-)) . (3.9) 
fc=i fc=i fc=i 

As a consequence, we can fix 03 small enough and L3 sufficiently large, so that 

4+i (*) - > ^(4+i - 4) > 

for any 1 < k < N — 1, when a < and L > L3. The first inequality in (|43p follows by writing 

a k +i(t) - a k {t) - a k +i(0) + a k (0) = / (a' k+1 (s) - a k (s)) ds > a*t. 

Jo 

Applying (|27p to the pair (rj(-,0),v(-,0)), we deduce the second inequality in (|4"3|) , 
Concerning ([44p . it follows from (|3 . 9|) that 



^2-K(t))2-^2-(c*)2| = 0(a) + 0((Lexp 



Therefore, one can decrease «3 and increase L3, if necessary, so that inequality (1441) holds for 
any 1 < k < N, when a < 03 and L > L3. This completes the proof of Corollary |H □ 

3.2 Monotonicity of localized scalar momentum 

This subsection is devoted to the almost monotonicity properties of suitably localized versions 
of the scalar momentum 

P(rj, v ) := - / rjv. 



2 . 

We derive them using the conservative form of the equation governing the integrand rjv of P. 
Lemma 3.1. Let (r],v) G C°([0, T], A/"V 3 (K)) be a solution to (IHUF1) . Then, 

dt[vv) = d x ((1 - 2 V )v 2 + ^ + (S 4 ( ^y ^ + + M 1 ~ *)) 011 [0,T]xR. (3.10) 

Proof. Under such regularity, it follows from (jHGPp that dtr] and c^v are continuous functions 
on [0, T] x M and that 

The conservation law given by (|3.10p then follows from the identity 

d" xV , {d x nf x 1 3 , _ /(3 - 2r / )(c\r ? ) 2 ■ 



^2(1 -r?) + 4(1- 77)2 ) - 2 ax[V + ln[L V)) °H 4(1 -77)2 J 



□ 



One can drop the smoothness assumptions in Lemma 13.11 by deriving an integral version of 
([3~TUD . 
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Corollary 3.1. Let (rj,v) E C°([0, T], MV(M)) be a solution to (iHCP]) . and /ei ip E C°([0,T], 
Cf(M))nC 1 ([0,T],C fe °(M)). T/ien, on [0,7] 

dt ( jf V^>) = J ftotov - Jjw ((i - 2^ 2 + ^ + (3 =^|^ )-\i d ^ + " ^) • 

77 (3.11) 

Proof. It follows by approximation of the initial datum using Theorem [TJ □ 

In particular, taking ijj = 1 in Corollary 13.11 yields the conservation of P. 
We turn now to localized versions of P. We deduce from Corollary 13.11 

Proposition 3.1. Let (rj,v) E C°([0, T],J\fV(R)) be a solution to dHGPl) . Lei ^ E Cf(R) 6e 
a non- decreasing function such that \d%SS?\ < Cod x ^ on R /or some positive number Cq. Let 
X E (^([OjT^M), t G [0,T] and r > be such that 2C r$ < I, and 



v(to,x) <- - An ^ _ 2 . , (3.12) 



1 *'(to) 

2 4(1 - 2C t 2 

/or any x E [X(io) — Ro,X(to) + Po] anc ^ some positive number Rq. Then, 
d t ( f Vfo(x-X(t)))ri(x,t)v(x,t)dx) <r (6 + — Wt/,v) sup d x *(x), (3.13) 

WM /|t=to t-r/ SU p/ |x|>r i?o 

w/iere 

r/sup := sup {?7(t, x), (i, x) E [0, T]xl}< 1. 

Proo/. Set J := [X(t ) - P , X(t ) + Po]- It follows from ([332"]) and the condition 2C t 2 < 1 
that T) < 1/2 on Lq, and real analysis then implies that 

|?7 + ln(l-77)| < 77 2 on J . (3.14) 

Applying (|3.11|) . we may rewrite the left-hand side of ()3.13p as 

- jf r ^*((l - 2^ 2 + ^ + X\to) V v + (3 ~^y 2 ) - i jf r 3 5^(ry + ln(l - r,)) , 

where <9 X \I' and d^> are evaluated at the point tq(x — X(to)). 

For x E Io> we deduce from (|3.14|) and the bound \d x ^\ < Cq8 x ^ that 

T d x V ((1 - 2 V )v 2 + ^ + X'(i V) + ^^/(t? + ln(l - r?)) 

> 7bfix*((l " 2r/)u 2 + (i - C7 t 2 )77 2 - \X'(t )\ V v) > 0, 

where the non-negativity of the quadratic form, which yields the last inequality, follows from 
assumption (|3.12p . Since r] < 1/2 on Lq, we also have 



i (3 -2r ] )(d x V ) 2 
4(1 -nf 



for x G Lq, and therefore, 

(3-2r ? )(d r?? ) 2 x , Tq 



r a^ ((1 - ^ + \+ X'(t )nv + ^ ) + ^fl£*fo + k(l "»?))< (3.15) 
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on Jo- 
lt remains to consider the case x ^ Jo. In view of the definition of r) snp , we have 

2 (1 -rj)v 2 , . o (1 - ri)v 2 

Tr^< (l-2^ 2 <4 l " , 



1 ??sup 2 

(3-2? ? )(a i .? ? ) 2 ^ r 6 "I (d x r]f 



- 4 1 - 77 2 " I ' 1-7WJ 8 1 



4(1 -r/) 2 " I 1 - r/ sup J 8(1 - rf) 
as well as 

|r? + ln(l - 77)1 < £ and \X' {t )r,v\ < 2^2^- + ^ ^ ~ ^ 



1 %up 4 4 1 ?7sup 

where we have used the inequality |X'(io)| < V2- Therefore, 

T d x V ((1 - 2 V )v 2 + ^ + X'(t )rtv) + \r^[ n + ln(l - rj)) 

<t (6 + - — )e(r],v) sup d x ^/(x) 

v J- — r/sup ' |a;|>roi?o 



(3.16) 



on 1R\ Iq, an d the conclusion follows from (|3. 15|) and (j3. 16)) by integration. □ 

We are finally in position to provide the 

Proof of Proposition^ Let to 6 [0, T\. We apply Proposition 13.11 with the choices to := ^0) 
ty(x) := (1 + (x))/2, ro := zv/16 and X(t) := (ak(t) + cik+±(t))/2. A quick computation first 
shows that for this choice of ty, the inequality \d x ^\ < Cq8 x ^ holds for Co := 4. The condition 
2CqTq < 1 is satisfied since v c * /16 < V2/16. We turn now to condition (|3.12p . It follows from 
(|44"|) . namely 

^-(a'^to)) 2 ^ and ^2-(4 +1 (t )) 2 >^, 

that 

X'(i ) 2 <2-^. (3.17) 

Using the inequality 1/(1 — s) < 1 + 2s for s 6 [0, 1/2], we infer from (|3.17p and the inequality 
v c */16 < 1/4 that 



2 4(1 - 8(v*) 2 ) ~ 2 4(1 - 8(z/*) 2 ) ~ 2 V2 16 / V v ; / ~ 32 ' 

At this stage, we come back to the decomposition 

(r/(-, t),v(-,t)) = R <t)Mt) (-) + e(-,t). 

In view of (|25|) . we have ||e(-,t)||x < K\a. Hence, by the Sobolev embedding theorem, we may 
choose 04 < as so that, if a < a^, then 

\H;t)U°° < (3-18) 
It follows from (|3.18p that (|3. 12[) holds everywhere on the set 

S := {x G R, s.t. |r7 c(to))a(to) (x)| < 
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Combining the definition of r/ c0 , the explicit form ([5]) of the solitons, condition (|43p on a(to) 
and the inequality z/ c > z^ c */2 stated in (|28p . it follows that the set S contains the interval 
[X(t )-R ,X(t ) + R } where 

(L-l) + a*t 2 /256AT 
i? := In I — — 

2 u c * \ z^r. 

Finally, since < d x ^(x) < 2exp(— 2|x|) for any x G R, we obtain from (|3,13p that 

l (Qi(t )) Mo£0 (e Xp (-^±^) 
and the proof is completed. □ 
We finally give the 

Proof of Corollary [JJ Since the energy E and the momentum P are conserved quantities for 
equation (jHGPp . we infer from definition (|45p that 



S (GW)=E(4-4«) S (%(*) 



fc=2 



Corollary [5] follows by combining Proposition [5] and the ordering condition (|42p . □ 

A A quantitative version of the implicit function theorem 

In this appendix, we give the proof of the following version of the implicit function theorem 
which we have invoked in order to establish Proposition |4j 

Proposition A.l. Let E and F be two Banach spaces, V, an open subset of F, andcp, a function 
of class C 2 from E x V to F . We assume that 

(i) there exist an open subset UofV and a map x of class C 1 from U to E such that 

4>(x(y),y) = 0, 

for any y G U . 

(ii) there exist a positive constant K, as well as a continuous isomorphism A y from F to F, 
such that 

d ?As(2/),y) = A y (ld + H y ), (A.l) 

and 

dx4>(x(y),y) = A y Ty, (A. 2) 

for any y G U . In (jA.ip . the operator norm of the continuous linear mapping H y is less 
than 1/2, whereas in (|A.2p . the operator norm ofT y is less than K. 

(Hi) the operator norm of d 2 (j)r Xj y-\ is uniformly bounded with respect to x G E and y G V. 

(iv) there exists a positive constant r such that 

B(y,r) C V, 

for any y G U . 
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Then, there exist two positive numbers 5 and A := AK + A, such that, given any y G U , there 
exists a map ir y G C 1 (B(x(y),5),V) such that, given any t G B(x(y),5), z = ir y (t) is the unique 
solution in B(y,A5) of the equation 

<f>(t,z)=0. 

Moreover, the map ir y is Lipschitz on B(x(y),5), with Lipschitz constant at most A. 

Remark A.l. When y is fixed in U, assumptions (i) and (ii) are enough to apply the implicit 
function theorem in order to construct a map ir y G C 1 (B(x(y), 5), V) as in Proposition lA.il The 
difficulty comes from the fact that, in this case, the constants 5 and K depend a priori on y. 

Proof. The proof is similar to the usual proof of the implicit function theorem. We fix y G U 
and we set 

§(x,z) := (<$>i(x,z),<f> 2 {x,z)) := (x,Ay 1 </)(x,z)-T y {x)). 

The function $ is well-defined and of class C 2 from E x V to E x F. Moreover, its differential 
at the point (x(y),y) may be written as 

'Id 



d ®(s(y),y) y Id + H{y) 

In view of assumption (ii), the operator norm of the maps d&(x(y),y) ~ Id 1S ^ ess tli an 1/2 for any 
y G U. Combining with assumption (Hi), we infer the existence of a positive number p < r such 
that the operator norm d§^ x ^ — Id is less than 3/4 for any (x, z) G B(x(y), p) x B(y, p). 

Given a point t G E, we next set 

V(x, z) := (t, -T y {t)) + (x, z) - <f>(x, z). 

The map ^ is well-defined on E x V, and Lipschitz on B(x(y),p) x B(y,p), with Lipschitz 
constant at most 3/4. Moreover, one can check that 

*(x(y),y) = {t,y-T y (t-x(y))), 

so that, letting A := 4 + AK, 

\\^{x(y),y) - (x(y),y)\\ ExF < |, 

when t G B(x(y), p/A). In this case, the map f is a contraction from B(x(y),p) x B(y,p) 
to B(x(y),p) x B(y,p). By the fixed point theorem, there exists a unique point (xt,Zt) in 
B(x(y),p) x B(y,p) such that ^(x t ,z t ) = (x*,^), i.e. 

xt = t and </>(i, z t ) = 0. 

At this stage, we set 5 := p/A < p < r, and 

7r y (i) := z t G V. 

In view of the definition of Zt, given two points t\ and t2 in B(x(y), 5), we have 

-%(^) = -r 2/ (tl-t 2 ) + (%(*l) - $2(*l,%0l) - 7Tj,(t 2 ) - $2(*2,7rj,(*2)), 

so that 

3 

- %(*2)||f < - t 2 \\ E + - t 2 ||E + hy{h) ~ TT y (t 2 )\\ F ). 

As a consequence, the map 7r y is Lipschitz on B(x(y),S), with Lipschitz constant at most AK+3 < 
A. The fact that ir y is of class C 1 on B(x(y),5) follows from the implicit function theorem. This 
completes the proof of Proposition IA.1I □ 
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